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O 

' A numerical analysis of several phases in the cold neutral two flavor dense quark matter 

bJQI is realized. In the analysis, besides the normal, single plane wave LOFF, and color super- 

^ , conducting 2SC/g2SC phases, we also include two gluonic phases (the gluonic cylindrical 

• phase II and the gluonic color-spin locked one). It is shown that these two gluonic phases 

are dynamically realized and, on the basis of this analysis, the phase diagram with respect 
00 ^ to the coupling strength in this medium is suggested. In particular, it is shown that the 

. gluonic phases are energetically favorable in a wide region of the parameter space and they 

can resolve the problem of the 2SC/g2SC chromomagnetic instability in the whole region 
where it exists. On the other hand, there also exists a window where the LOFF phase is 
stable. The 2SC state and the normal one are realized in the strong and weak coupling 
regimes, respectively, as was expected. 

^: 
I) : 

^ , §1. Introduction 



■ In this paper, we will present a self-consistent numerical analysis of several phases 
\ in the cold neutral two flavor dense quark matter. In particular, we analyze gluonic 
I phases introduced in Ref. 1) and described recently in detail in Ref. 2). On the basis 

■ of this analysis, it is shown that a) certain gluonic phases connected with a first order 
\ phase transition are dynamically realized and b) the phase diagram with respect to 
' the coupling strength in this medium is suggested. In particular, it is shown that the 

gluonic phases are energetically favorable in a wide region of the parameter space 
i and they can resolve the problem of the 2SC/g2SC chromomagnetic instability in 

the whole region where it exists. 

It is expected that at sufficiently high baryon density, cold quark matter should 
\ be in a color superconducting state (for reviews, see Ref. 3)). On the other hand, 

' it was suggested long ago that quark matter might exist inside the central region 

of compact stars. This is one of the main reasons why the dynamics of the color 
superconductivity has been intensively studied. 

Bulk matter in compact stars should be in /3-equilibrium, providing by weak 
interactions, and be electrically and color neutral. The electric and color neutrality 
conditions play a crucial role in the dynamics of quark pairing. ^^^^^^ Also, in the 
dense quark matter, the strange quark mass cannot be neglected. These factors lead 
to a mismatch 5/i between the Fermi momenta of the pairing quarks. 

The mismatch 5^ tends to destroy the gapped (2SC) and gapless (g2SC) two- 
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flavor color superconducting phases and the color-flavor locked (CFL) and gapless 
color-flavor locked (gCFL) three-flavor phases. Establishing the genuine ground state 
(vacuum) for realistic values of 6fi in dense QCD is one of the central problems in the 
field. Besides the gluonic phases/)'^-* a number of other candidates for the ground 
state have been proposed. ^^^^^^^ 

An important factor in solving this problem is removing instabilities observed 
in gluon (color plasmon) channels in the 2SC/g2SC and gCFL phases. ^^^"^^^ As was 
revealed in Ref. 20), the 2SC/g2SC phases suffer from a chromomagnetic instability 
connected with the presence of imaginary Meissner masses of gluons. While the 8th 
gluon has an imaginary Meissner mass only in the g2SC phase, with the diquark 
gap A < 5fi (an intermediate coupling regime), the chromomagnetic instability for 
the 4- 7th gluons appears also in a strong coupling regime in the 2SC phase, with 
6ijl < A < \f2b[i. Later a chromomagnetic instability was also found in the three- 
flavor gCFL phase. 

Meissner and Debye masses are screening (and not pole) ones. It has been re- 
cently revealed in Ref. 22) that the chromomagnetic instabilities in the 4- 7th and 
8th gluonic channels in the 2SC/g2SC phases correspond to two very different tachy- 
onic spectra of plasmons. It is noticeable that while (unlike the Meissner mass) the 
(screening) Debye mass for an electric mode remains real for all values of b\i both in 
the 2SC and g2SC phases, ^'^^ the tachyonic plasmons occur both for the magnetic 
and electric modes. ■^^^ The latter is important since it clearly shows that this in- 
stability is connected with vector-like excitations: Recall that two magnetic modes 
correspond to two transverse components of a plasmon, and one electric mode cor- 
responds to its longitudinal component. This form of the plasmon spectrum leads 
to the unequivocal conclusion about the existence of vector condensates of gluons in 
the ground state of two flavor quark matter with A < \/^6fJ,, thus supporting the 
scenario with gluon condensates (gluonic phases) proposed in Refs. 1),2). 

Besides the instabilities in gluon channels, there exists an instability reflected in 
a negative velocity squared < of the physical diquark excitation (diquark Higgs 
mode) in the g2SC phase'^^^ (a similar instability was discussed in Refs. 24), 25)). 
This problem should be also resolved in the genuine ground state. 

In Refs. 1),2), by using the Ginzburg-Landau (GL) approach in the hard dense 
loop (HDL) approximation, it was demonstrated the existence of a particular gluonic 
phase in the vicinity of the critical point 6fi ~ A/\/2. Since the GL approach works 
only for second order phase transitions, a specific ansatz, leading to such a transition, 
was considered. It is: 

/.8^^5(4), f^3^9{Al), B^g{Al), C^g{Al), (M) 

where ^^'s are gluon fields and g denotes the QCD coupling constant. In this 
phase, the symmetry SU{2)c x [/(l)em x 5'0(3)rot is spontaneously broken down to 
50(2)i.ot- Here SU{2)c and U{l)em are the color and electromagnetic symmetries 
in the 2SC state, and SO {3)rot is the rotational group (recall that in the 2SC state 
the color SU{3)c is broken down to SU{2)c). In other words, the solution (jl-ip 
describes an anisotropic medium in which the color and electric superconductivities 
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coexist. It is noticeable that there are chromoelectric field strength condensates 
{E1) and (-EJ) in this phase, reflecting its non-Abelian nature. Since solution (jl-ip 
is cylindrically symmetric, it was called the gluonic cylindrical phase I in 2) (for the 
gluonic cylindrical phase II, see below). 

The first important step in realizing a numerical analysis in the gluonic phase 
outside the scaling region around 5/i ~ Ajypl was done in Refs. 26), 27). In that 
analysis, some approximations were made: the color chemical potentials /is and /ig 
were ignored and it was assumed that the mismatch bii and the diquark gap A in 
the gluonic phase are the same as those in the 2SC/g2SC one. However, it is unclear 
whether these approximations are justified in the gluonic phase. For example, while 
in the 2SC/g2SC phases /xs = and /xg is small, it is questionable whether they 
can be neglected in the gluonic one. As to (5/i and A in the gluonic phase, nothing 
is known about them away from the critical point either. 

The main goal of the present paper is a numerical self-consistent analysis of the 
gluonic phases. Unfortunately, at present it is very hard to realize such a study in the 
phase (jl-ip . The point is that as was shown in Ref. 2), outside the scaling region, two 
additional gluon condensates, /lx4 (= g{A^) and /is (= g^J^Qj^^ cannot be ignored. 
With the four condensates in Eq. (jl-ip and these two ones, plus two equations for 
b\i and A^ the numerical analysis in this phase is very involved. Fortunately, as was 
pointed out in 2) (based on the description of possible symmetry breaking structures 
in dense QCD), there may exist other gluonic phases (presumably connected with 
first order phase transitions). In some of them, it is easier to realize a numerical 
analysis than in the gluonic cylindrical phase I. 

Here we will be interested in the following two ones. First, a phase with C = 
in Eq. (jl-ip is not excluded: 

/X8 = ^5(^8), /.3=5(^g), B = g^Al). (1-2) 

It was called the gluonic cylindrical phase 11.^^ In this phase, while the rotational 
50(3) is also spontaneously broken down to SO(2)^o\.-, the electromagnetic ^7(1) is 
preserved. We emphasize that since a chromoelectric field strength condensate {E\) 
is generated in this phase, it also possesses a non-Abelian nature. 

The second phase that will be analyzed here is the so called gluonic color-spin 
locked (GCSL) phase with the condensates 

^8 = ^5(4), K^g{A\^=g{A\). (1-3) 

In this phase, the symmetry SU{2)c x U{l)em x SO{3)rot is spontaneously broken 
down to the SO {2)dis,g one, with the generator which is an appropriate linear com- 
bination of C7^ in the color SU (2)^ and the generator of spatial y — z rotations. It is 
noticeable that there exist both chromoelectric and chromomagnetic field strengths 
condensates in the GCSL phase. Also, since the electromagnetic U{1) is sponta- 
neously broken in this phase, it (like the gluonic cylindrical phase I) describes an 
anisotropic superconductor. 
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It is quite nontrivial whether or not the gluonic phases (ll-2p and (I1-3P connected 
with the first order phase transition are dynamically reahzed. Fortunately, there 
are not so many parameters in these phases that will allow us to perform numerical 
calculations for solving self-consistently both the neutrality conditions for fj.^, ^ug, 5/i 
and the gap equations for A, B and/or K (the five coupled equations in the gluonic 
cylindrical phase II and the four coupled equations in the GCSL phase). 

As a benchmark, we will also study the single plane LOFF phase, ^^^"^^-^ which 
is gauge equivalent to^^'^°^ 

H, = ^9{Al), q^-l^g{Al). (1-4) 

We will determine which of the phases 2SC/g2SC and ([F2]), dliS]), and ([HI) is 
energetically favorable in the bulk of the parameter space. 

It turns out that there indeed exist dynamical solutions corresponding to the 
gluonic cylindrical phase II and the GCSL one. Actually, these phases are ener- 
getically more favorable than the single plane wave LOFF and 2SC/g2SC phases 
in a wide region of the parameter space. For realistic values // = 400 MeV and 
A = 635.5 MeV, where /i and A are the quark chemical potential and the cutoff in 
the (gauged) Nambu-Jona-Lasinio (NJL) model, the gluonic cylindrical II and GCSL 
phases are more stable than the normal, LOFF and g2SC/2SC phases in the wide 
interval 6.7 x 10^ MeV < Aq < 1.6 x 10^ MeV. Here we introduced the 2SC gap 
parameter Aq defined at (5/U = (it essentially corresponds to the strength of the 
four-diquark-coupling constant in the (gauged) NJL model). 

On the other hand, the neutral normal phase with no condensates is the ground 
state at Aq < 6.5 x 10^ MeV (weak coupling regime) while the single plane LOFF 
phase is favorable in the narrow window 6.5 x 10^ MeV < Z\o < 6.7 x lO^MeV. At 
^0 > 1.6 X 10^ MeV (strong coupling regime), the 2SC phase is realized. 

A noticeable point is that the gapless g2SC solution exists only at 92.2 MeV < 
^0 < 134.6 MeV, i.e., it occurs later (at larger values of Aq) than the gluonic 
phases. This strongly suggests that the gluonic phases can resolve the problem of 
the chromomagnetic and plasmon instabilities in the 2SC/g2SC phases. 

§2. Model 

As in Refs. 1), 2), we will use the gauged Nambu-Jona-Lasinio (NJL) model 
with two light (up and down) quarks, whose current masses will be neglected. The 
Lagrangian density is given by 

C = ^Piip + /2o7°)V' + Ga 
where 

D, ^d,- igA'lT'^, f;, ^ d^Al - d^A^ + gf^^A'^Al (2-2) 

with £ = e*-' and e" = e"^''' are the totally antisymmetric tensors in the fiavor and 
color spaces, respectively; i,j = u,d and a,/3, 7 = r (red), g (green), b (blue). 



(V;^iee°75^)(?/^iee°75V''^') 



(2-1) 
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The matrices T^'s are the generators of SU{2>) in the fundamental representation, 
and f°-^'^^s are the structure constants of SU{3). The charge-conjugate spinor V'*^ 
is ip^ = C'ljF with C = ^7^7*^. We do not introduce photon field and therefore 
the electric charge is connected with a global U{\) symmetry in the model. On the 
other hand, the whole theory contains free electrons, although we do not show them 
explicitly in Eq. ()2-ip . In /3-equilibrium, the chemical potential matrix connected 
with the conserved global quark and electric charges, is 

/^o = - ^eQcm = Ao - (^^T-a- (2-3) 

Here /U and are the quark and electron chemical potentials, /io = /i — ^ and 
6fi = ^ (the baryon chemical potential fiB is IJ-b = 3//.) The matrix ra = diag(l, —1) 
acts in the flavor space. Also, in the 2SC/g2SC phases, the color neutrality condition 
requires a non-zero color chemical potential /^s-^^^ 

The following comment is in order. In the full theory, the four-quark coupling 
constant Ga (or Aq) is not a free parameter but expressed through the QCD coupling 
g and the quark chemical potential /i. Needless to say, that at present this relation 
for ^0 is unknown. Because of that, we will treat as a free parameter, allowing it 
to vary in a wide interval that includes its typical values considered in the literature. 

Let us introduce the diquark field ~ i-ip'-' ee'^j^ip. We will use a gauge in 
which 

= 0, = 0, = A, (2-4) 

where the anti-blue component A is real. The gap A is given by the vacuum expec- 
tation value (VEV) of Z\, ^ = {A). 

In the presence of both diquark field A and gluon fields by using the Nambu- 
Gor'kov spinor 

^.(/c). (2.5) 

the inverse fermion propagator S^^ can be written as 

SiHP) = {^^' (2-6) 

with 

KgrHP) ^ (po + /^o)7°-7•p + 54"T^ (2-7) 

[GoJ'HP) ^ (po-/^o)7°-7-p-54'^(rT, (2-8) 

A- = -ieeSs^, ^+ = 7°(/^^)S° = -ieeSs^, (2-9) 

where = {pojP) is the energy-momentum four vector. Note that in this case the 
color neutrality requires in principle VEVs of the time components of all gluon fields: 

^i,2,...,7 = 5(^S'''-''), l^8^^9{Al). (2-10) 
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We will utilize the fermion one-loop approximation, closely related to the HDL 
one, in which only the dominant one-loop quark contribution is taken into account 
while the contribution of gluon loops is neglected. On the other hand, we keep the 
tree contribution of gluons. This is because we want to check how sensitive are the 
results to the choice of the value of the QCD coupling g. 

In this approximation, the effective potential including both gluon and diquark 
fields is: 

(2.11) 

where we also added the free electron contribution. The following subtraction scheme 
will be used: 

+1 1 ^^TrlnS-\P;fx = Sfi = iXa = A = 0,A''), (2-12) 
in which a spurious mass term yl^lA^p is removed. 

§3. Numerical analysis of the effective potential with gluon condensates 

3.1. Method 

We will perform the numerical analysis by following the approach developed in 
Ref. 7). The essence of this approach is the following. In order to calculate the 

fermion determinant, it is sufficient to know the energy eigenvalues of the inverse 
fermion propagator. Then, after integrating over pQ in the loop integral, we obtain 



a=l 



(3-1) 



where EaS are the energy eigenvalues of a n x n matrix S~^{P) including both 
positive (particle) and negative (anti-particle) eigenvalues. 

Before explicit calculations, let us show how one can block-diagonalize the inverse 
propagator in the flavor and chirality spaces: this will allow to reduce the 48 x 48 

matrix to a 12 x 12 one. We write: 



with 



det Sg^ = det 



det 



e 
1 



Sg ^ 



1 

e 



T -> 
■1, 



(3-2) 
(3-3) 

(3-4) 
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Fig. 1. The diquark gap A for the gluonic 
cylindrical phase II and the 2SC/g2SC 
phases. At Zio = 134.6 MeV, the g2SC 
phase turns into the 2SC one. Here the 
values A = 653.3 MeV and = 400 MeV 
were used. 
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Fig. 2. The mismatch &\i in the gluonic cylin- 
drical phase II and the 2SC/g2SC phases. 
At ^0 = 134.6 MeV, the g2SC phase turns 
into the 2SC one. The values A — 653.3 
MeV and n = 400 MeV were used. 



where we defined 

A^eoi = A = gA''T'', fl^ = flo + 5fiT3. (3-5) 

The expression (|3-3p is flavor-diagonal. 

Now, since the current quark masses are neglected, the theory is chiral invariant. 
Therefore one can decompose the inverse propagator into the right- and left-handed 
parts. This reduces the determinant of the 48 x 48 matrix Sg^ to that of two 12 x 12 
flavor-diagonal matrices. For the single plane wave LOFF phase, we can further 
reduce the size of the matrix to a 4 x 4 one. 

3.2. Results 

The form and the structure of the the effective potential (l2T2p were analyzed 
numerically. We search for stationary points of the potential, which correspond to 
the solutions of the neutrality conditions and the gap equations. In the analysis, we 
take realistic values // = 400 MeV and A = 653.3 MeV. The results are qualitatively 
unchanged in the region of the parameters with 300 MeV < fj, < 500 MeV and 
653.3 MeV < yl < 1 GeV. 

We find that the gluonic cylindrical phase II does exist in the region 

6.5 X 10^ MeV < < 1-6 x 10^ MeV . (3-6) 

In this case, the tree gluon term in Eq. ()2-lip is negligible, so that the results are 
not sensitive to the choice of the value of the strong coupling as = g'^ /A-k. We depict 
the solutions for A and 5^, following from the neutrality conditions, in Figs. [T]and 
[21 respectively. Note that due to the neutrality conditions, the values of A and 5^ in 
the gluonic cylindrical phase II is significantly different from those in the 2SC/g2SC 
ones. The behavior of the gluon condensate B is presented in Fig. O The phase 
transition from the gluonic cylindrical phase II to the 2SC one seems to be a weakly 
first order one. This conclusion is, however, not definite, because the accuracy is 
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not so good near this critical point. On 
the other hand, there definitely exists a 
first order phase transition between the 
cylindrical phase and the normal one: 
as one can see in Fig. [3l the value of 
the condensate B is large at the starting 
point of the former. 

In Fig. m the neutral effective po- 
tential at ^0 = 110 MeV is shown, with 
5fj,, /U3 and fig determined from the neu- 
trality conditions. It is clear that the so- 
lution for the gluonic cylindrical phase 
II corresponds to a minimum of the neu- 
tral effective potential. 

The values of the color chemical po- 
tentials ^3 and /is is found to be at most l/^sl ~ 0(50 MeV) and l/^sl ~ 0(40 MeV) 
for 300 MeV < < 500 MeV and 653.3 MeV < yl < 1 GeV. At the starting point 
of this phase, the value of the gap A is zero. At its endpoint, the ratio of 5fi/A is 
6i^/A = 0.84-0.86 > 1/V2 = 0.707 for these values /i and A. This result is consistent 
with that in Ref. 32). 

Let us now turn to the results of the analysis for the GCSL phase. It exists in 
the region 

7.6 X 10^ MeV < < 1-6 x 10^ MeV, (3-7) 

for as = 1. Note that the GCSL phase starts at a larger than the gluonic 
cylindrical phase II. The values of A and 6fj, in these two phases are not so different. 
However, unlike the case of the gluonic cylindrical phase II, the gap A at the starting 
point in the GCSL phase is quite big, which is a strong indication on a first order 
phase transition between this phase and the normal one. 

The form of the solution for the gluon condensate K is depicted in Fig. [3l Note 
that there is a tiny discrepancy (~ 0.5%) between the values of the endpoints of 
these two gluonic phases. A more accurate analysis is required to answer whether it 
is a real thing or an artifact. On the other hand, a large value of the condensate K 
at the starting point strongly supports the conclusion made above that a first order 
phase transition takes place between the GCSL phase and the normal one. 

Let us now compare the free energy for the 2SC/g2SC, LOFF, GCSL and glu- 
onic cylindrical phases. The corresponding numerical results are shown in Fig. [5j In 
the weak coupling regime with Aq < 65 MeV, the normal phase with no condensates 
is realized. The single plane wave LOFF phase is the most stable one in a narrow 
window with 65 MeV < Aq < 67 MeV. The gluonic cylindrical/GCSL phases are 
energetically favorable in a wide region with 67 MeV < Aq < 1.6 x 10^ MeV, corre- 
sponding to the intermediate and strong coupling regimes. At Aq > 1.6 x 10^ MeV 
the 2SC phase is realized. 

Comments concerning the sensitivity of the dynamics in the GCSL phase with 
respect to the strong coupling Og are now in order. In the region close to the endpoint 
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Fig. 3. The gluon condensates B and K in the 
gluonic cylindrical phase II and the GCSL 
phase. The values A = 653.3 MeV, fj. = 400 
MeV, and Us = 1 were used. 
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of this phase, its free energy V^s is es- 
sentiaUy independent of a^, i.e., the tree 
gluon term is irrelevant there. On the 
other hand, the value of at the start- 
ing point is rather sensitive to the choice 
of Us- For example, 

Z'o*'''*(GCSL) = 82, 76, 69 MeV (3-8) 

for Us = 0.75, 1.00, 1.25. The sensitiv- 
ity of the free energy V^g to the choice 
of as is also higher in that region, al- 
though it is not too strong: the varia- 
tions of V^ff are at most 0.4 MeV/fm^ for 
Us = 0.75-1.25. The latter is compara- 
ble with the difference between the free 
energies of the LOFF and GCSL phases 
in the region Aq < 95MeV at = 1 
(see Fig. [5]). The origin of the sensi- 
tivity in that region is connected with 
a contribution of the tree gluonic term 
in the effective potential (|2-lip . which 
is provided by a chromomagnetic field 
strength condensate. At the starting 
point, it is g{Fy^) = ~ (110-130 MeV)^ for = 0.75-1.25. It is impor- 

tant that, despite this sensitivity, qualitative results of the analysis in this region, 
such as a strong first order phase transition at the starting point, etc., are robust. 

As one can see in Fig. [SJ there exists a critical value Aq in the region 76 MeV < 
^0 < 1-6 X 10^ MeV dividing the gluonic cylindrical and GCSL phases: the latter 
becomes more stable at > ^o"^- However, since the value of ^q"^ is sensitive to 
the choice of a^, we do not give here a concrete number for it. 

It is instructive to compare our results with those in Ref. 27), where (in the 
present language) the gluonic cylindrical phase II was also numerically analyzed, 
but by utilizing the approximations described in the Introduction. In both these 
analyses, this gluonic phase is energetically more favored than the single plane wave 
LOFF one in a wide range of the parameter space. On the other hand, unlike 27), 
the gluonic phase in our analysis appears essentially earlier than the g2SC phase 
which exists only in the region 92.2 MeV < Aq < 134.6 MeV. This feature is 
qualitatively unchanged for 300 MeV < < 500 MeV and 653.3 MeV < yl < 1 GeV. 
It strongly suggests that the gluonic phase can resolve the problem of the 2SC/g2SC 
chromomagnetic instability in the whole region where it exists. 

§4. Summary and discussions 

The main result of this paper is establishing the existence of the gluonic cylin- 
drical phase II and the GCSL one that led us to the suggestion of the form of the 
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Fig. 4. The contour plot of the neutral effec- 
tive potential for the gluonic cylindrical 
phase II. Since the neutrality conditions for 
Sfi, ^3 and fis were solved at each point 
{A,B), their values are different at differ- 
ent points. The reference point for the free 
energy is its value in the normal phase. The 
values Ao = 110 MeV, A = 653.3 MeV, and 
^ = 400 MeV were used. 
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phase diagram in the cold two flavor 
dense quark matter. These phases are 
energetically more favorable than the 
2SC/g2SC and single plane wave LOFF 
ones in a wide region of the parameter 
space. This indicates that the gluonic 
phases can resolve the problem of the 
chromomagnetic and plasmon instabili- 
ties in two-flavor dense QCD. We also 
found that the single plane wave LOFF 
state is energetically favorable in a win- 
dow 65 MeV < ^0 < 67 MeV. In this 
region, the LOFF state is probably free 
from the chromomagnetic instability.^'^) 
It would be also worthwhile to examine 
the multiple plane wave LOFF state. ^^-^ 
The 2SC state and the normal one are 
realized in the strong and weak coupling 
regimes, respectively, as was expected. 

There is still quite a number of 
questions that need to be answered. It 
is still a nontrivial problem whether or not the < instability in the diquark 
Higgs mode^^)"^^) is absent in the gluonic phases. It would be also important to 
calculate explicitly the Meissner masses squared in order to prove the absence of 
the chromomagnetic instability there. Note that their values are essentially differ- 
ent from the eigenvalues of the curvature of the neutral effective potential. Indeed, 
while the former are calculated at fixed values of the chemical potentials 5^, /X3 and 
/xg (corresponding to the given gap A and the condensates B or K), the latter are 
calculated by varying the effective potential when these chemical potentials are also 
varying with A and B (or K) through the neutral conditions. This problem will be 
considered elsewhere. 

The dynamics with vector condensates of gauge fields, such as in the gluonic 
phases, are very rich and one can expect intriguing phenomena there, such as the 
existence of exotic hadrons,^-*'^) roton-like excitations,^^) and vortex-like ones.^^^ 
The GCSL phase is especially interesting: it describes an anisotropic superconductor. 
Moreover, we expect an abnormal number of gapless Nambu-Goldstone (NG) bosons 
in its spectrum. This expectation is based on recent results in Ref. 35) obtained in 
the gauged a-model with chemical potential for hypercharge. In that work, a phase 
with vector condensates, whose properties are very similar to those in the GCSL one, 
has been revealed (it was called the gauge-flavor-spin locked phase). It was shown 
that there is an abnormal number of NG bosons in that phase. We hope to turn 
to these dynamical issues elsewhere. 



*' It is the same phenomenon as that found in a non-gauge relativistic model at finite density 
in Ref. 36). 
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Ao (MeV) 

Fig. 5. The free energy in different phases. 
Its value in the normal phase is chosen as 
the reference point. The bold and dashed 
curves correspond to the gluonic cylindri- 
cal and GCSL phases, respectively. The 
dotted and dash-dotted ones correspond to 
the single plane wave LOFF phase and the 
2SC/g2SC ones. The values A = 653.3 
MeV, = 400 MeV, and = 1 were used. 
In the inset, the free energy in the region 
near the normal phase is shown. 



Gluonic phases and phase diagram in neutral two flavor dense QCD 



11 



Acknowledgments 

M.H. acknowledges useful discussions with H. Abuki. A part of this work was 
done when VA.M. visited Yukawa Institute for Theoretical Physics, Kyoto Univer- 
sity. He is grateful to Professor Taichiro Kugo and Professor Teiji Kunihiro for their 
warm hospitality. The research of M.H. was supported by the JSPS Grant-in-Aid 
for Scientific Research (B) 18340059. The work of V.A.M. was supported by the 
Natural Sciences and Engineering Research Council of Canada. 

References 

1) E. V. Gorbar, M. Hashimoto and V. A. Miransky, Phys. Lett. B 632 (2006), 305. 

2) E. V. Gorbar, M. Hashimoto and V. A. Miransky, Phys. Rev. D 75 (2007), 085012. 

3) K. Rajagopal and F. Wilczek, in Handbook of QCD, M.Shifman, ed., (World Scientific, 
Singapore, 2001), hep-ph/0011333'. 

M. G. Alford, Ann. Rev. Nucl. Part. Sci. 51 (2001), 131. 
D. K. Hong, Acta Phys. Polon. B 32 (2001), 1253. 
S. Reddy, Acta Phys. Polon. B 33 (2002), 4101. 
T. Schafer, hep-ph/0304281 

D. H. Rischke, Prog. Part. Nucl. Phys. 52 (2004), 197. 

M. Buballa, Phys. Rep. 407 (2005),' 205. 

M. Huang, Int. J. Mod. Phys. E 14 (2005), 675. 

I. A. Shovkovy, Found. Phys. 35 (2005), 1309. 

4) D. Ivanenko and D. F. Kurdgelaidze, Astrofiz. 1 (1965), 479; Lett. Nuovo Cim. 2 (1969), 
13. 

N. Itoh, Prog. Theor. Phys. 44 (1970), 291. 

F. lachello, W. D. Langer and A. Lande, Nucl. Phys. A 219 (1974), 612. 
J. C. Collins and M. J. Perry, Phys. Rev. Lett. 34 (1975), 1353. 

5) K. lida and G. Baym, Phys. Rev. D 63 (2001), 074018 [Errata; D 66 (2002), 059903]. 

6) M. Alford and K. Rajagopal, JHEP 0206 (2002), 031. 

7) A. W. Steiner, S. Reddy and M. Prakash, Phys. Rev. D 66 (2002), 094007. 

8) M. Huang, P. f. Zhuang and W. q. Chao, Phys. Rev. D 67 (2003), 065015. 

9) L Shovkovy and M. Huang, Phys. Lett. B 564 (2003), 205. 
M. Huang and L Shovkovy, Nucl. Phys. A 729 (2003), 835. 

10) H. Abuki, M. Kitazawa and T. Kunihiro, Phys. Lett. B 615 (2005), 102. 
H. Abuki and T. Kunihiro, Nucl. Phys. A 768 (2006), 118. 

11) S. B. Ruster, V. Werth, M. Buballa, I. A. Shovkovy, and D. H. Rischke, Phys. Rev. D 72 
(2005), 034004. 

12) D. Blaschke, S. Fredriksson, H. Grigorian, A. M. Oztas and F. Sandin, Phys. Rev. D 72 
(2005), 065020. 

13) M. G. Alford, J. A. Bowers and K. Rajagopal, Phys. Rev. D 63 (2001), 074016. 

For a review, see, e.g., R. Casalbuoni and G. NarduUi, Rev. Mod. Phys. 76 (2004), 263. 

14) J. A. Bowers and K. Rajagopal, Phys. Rev. D 66 (2002), 06;5002. 
M. MannareUi, K. Rajagopal and R. Sharma, hep-ph/070 20'2T] 

15) S. Reddy and G. Rupak, Phys. Rev. C 71 (2005), 025201. 

16) M. Huang, Phys. Rev. D 73 (2006), 045007. 

17) A. Kryjevski, hep-ph/0508180 

T. Schafer, Phys. Rev. Lett. 96 (2006), 012305. 

A. Gerhold, T. Schafer and A. Kryjevski, Phys. Rev. D 75 (2007), 054012. 

18) X. B. Zhang and J. I. Kapusta, Phys. Rev. D 75 (2007), 054028. 

19) R. Gatto and M. Ruggieri, Phys. Rev. D 75 (2007), 114004. 

20) M. Huang and L A. Shovkovy, Phys. Rev. D 70 (2004), 051501(R); ibid. D 70 (2004), 
094030. 

21) R. Casalbuoni, R. Gatto, M. MannareUi, G. Nardulh and M. Ruggieri, Phys. Lett. B 605 
(2005), 362 [Errata; B 615 (2005), 297]. 

M. Alford and Q. Wang, J. Phys. G 31 (2005), 719. 



12 



M. Hashimoto, V. A. Miransky 



K. Fukushima, Phys. Rev. D 72 (2005), 074002. 

22) E. V. Gorbar, M. Hashimoto, V. A. Miransky and I. A. Shovkovy, Phys. Rev. D 73 (2006), 
111502(R). 

23) M. Hashimoto, Phys. Lett. B 642 (2006), 93. 

24) K. lida and K. Fukushima, Phys. Rev. D 74 (2006), 074020. 

25) I. Giannakis, D. Hou, M. Huang and H. c. Ren, Phys. Rev. D 75 (2007), 011501(R); ibid. 
D 75 (2007), 014015. 

26) K. Fukushima, Phys. Rev. D 73 (2006), 094016. 

27) O. Kiriyama, D. H. Rischkc and I. A. Shovkovy, Phys. Lett. B 643 (2006), 331. 

28) A. Gerhold and A. Rebhan, Phys. Rev. D 68 (2003), 011502(R). 

D. D. Dietrich and D. H. Rischke, Prog. Part. NucL Phys. 53 (2004), 305. 

29) L Giannakis and H. C. Ren, Phys. Lett. B 611 (2005), 137; NucL Phys. B 723 (2005), 
255. 

L Giannakis, D. f. Hou and H. C. Ren, Phys. Lett. B 631 (2005), 16. 

30) E. V. Gorbar, M. Hashimoto and V. A. Miransky, Phys. Rev. Lett. 96 (2006), 022005. 

31) L. He, M. Jin and P. Zhuang, Phys. Rev. D 75 (2007), 036003. 

32) O. Kiriyama, Phys. Rev. D 74 (2006), 074019; ibid. D 74 (2006), 114011. 

33) V. P. Gusynin, V. A. Miransky and L A. Shovkovy, Phys. Lett. B 581 (2004), 82; Mod. 
Phys. Lett. A 19 (2004), 1341. 

34) E. V. Gorbar, J. Jia and V. A. Miransky, Phys. Rev. D 73 (2006), 045001. 

35) A. Buchel, J. Jia and V. A. Miransky, Nucl. Phys. B 772 [FS] (2007), 323. 

36) V. A. Miransky and I. A. Shovkovy, Phys. Rev. Lett. 88 (2002), 111601. 

T. Schafor, D. T. Son, M. A. Stephanov, D. Toublan and J. J. M. Verbaarschot, Phys. 
Lett. B 522 (2001), 67. 



